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ABSTRACT 

4U 2206+54 is a high mass X-ray binary whic h has been suspect ed to contain a neutron star accreting 
from the wind of its companion BD +53° 2790. iReig et all ((20091 ) have recently detected 5560 s period 
pulsations in both RXTE and INTEGRAL observations which they conclude are due to the spin of 
the neutron star. We present observations made with Suzaku which are contemporaneous with their 
RXTE observation of this source. We find strong pulsations at a period of 5554 ± 9 s in agreement 
with their results. We also present a reanalysis of BeppoSAX observations of 4U 2206+54 made in 
1998, in which we find strong pulsations at a period of 5420 + 28 seconds, revealing a spin-down trend 
in this long-period accreting pulsar. Analysis of these data suggests that the neutron star in this 
system is an accretion-powered magnetar. 

Subject headings: accretion, accretion disks — X-rays: binaries — pulsars: individual (4U 2206+54) 
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1. INTRODUCTION 

First observed by Uhuru and Ariel V, the low luminos- 
ity (~10 35 erg cm~ 2 s _1 ) hard x-ray source 4U 2206+54 is 
one of a handful of well studied galactic high mass x-ray 
binaries for which the nature of the accreting compact 
star is uncertain. 

Optical and UV spectroscopy have showed the optical 
counterpart, BD +53° 2790 to be an unusual 09 active 
sta r with strong wind reson ance lines in the ultravio- 
let (|Negueruela &; R cig 2001ft . Several authors have con- 
cluded from the high column depth and the flaring ob- 
served in the x-rays that the compact object is accreting 
from a strong stellar wind. Ribo et al. (2006) determine 
from the an IUE spectra a relatively low wind velocity 
of -350 km s _1 . 

Analysis of the first five years of the RXTE ASM light- 
curve of 4U 2206+54 revealed a periodic modulation at a 
period of 9.57 days, which w as interpreted as the o rbital 
period of the binary system (jCorbet &: P eclc 2001). Ra- 
dial velocity measurements were made for BD +53° 2790, 
however it was not possible t o independ ently determine 
the orbital period from them (Blav 2005) . Recent analy- 
ses of the Swift /BAT and RXTE/ A SM light-curves hav e 
thrown this picture into confusion (jCorbet et al.ll2007t ). 
The recent data shows a modulation with a period of 
19.25 days - twice what was previously thought to be 
the period. If 19.25 days is the orbital period, then the 
profile of the orbital modulation has evolved from having 
two peaks per orbit to one peak per orbit. 

Using observations extended over 7 days from the Pro- 
portional Counter Ar ray (PCA)on the Rossi X-ray Tim- 
ing Explorer (RXTE) IReig et al.l(|2009f ) discovered pulsa- 
tions at a period of 5559+3 s. Detecting pulsations in this 
period range from a low Earth orbit is usually difficult be- 
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cause the pulse period is near the orbital period, resulting 
in poor sampling of pulses. Their data however include 
several intervals which were uninterrupted by Earth oc- 
cultation. They also confirmed this discovery using ob- 
servations from IBIS on the International Gamma-ray 
Astrophysics Laboratory (INTEGRAL) which has a high 
Earth orbit with a 3 day period. 

In section 2 we present results from new observations 
of 4U 2206+54 we have made with Suzaku. These ob- 
servations were c ontemporaneous with those of RXTE 
(|Reig et al.ll2009ft . We find strong pulsations at a pulse 
period of 5554+9 s, in agreement with the RXTE results. 
In section 3 we present a reanalysis of the BeppoSAX ob- 
servations of 4U 2206+54 which show strong pulsations 
at a period of 5420 ± 28 s, in section 4 we discuss simu- 
lations of Suzaku and BeppoSAX data, in section 5 the 
EXOSAT observations of 4U 2206+54, and in section 6 
we briefly analyze the evolutionary status of the system 
and the nature of its primary component. We conclude 
that the primary component of the system is a neutron 
star, which has a huge (~ (3 — 5) x 10 15 G) surface field 
and is accreting material onto its surface from a disk. In 
this light 4U 2206+54 appears to be the first accretion- 
powered magnetar identified so far. 

2. ANALYSIS OF THE SUZAKU LIGHT-CURVE 

We obtained an observation of 4U 2206+54 with 
Suzaku with the primary objective of searching for long- 
perio d pulsations with the X-ray Imaging Spectrometer 
(XIS; iKomava et al.ll2007l ) . To minimize the number and 
size of data gaps the observation was time constrained to 
a phase of the space-craft precession cycle were the source 
was always above the Earth horizon. The observation oc- 
curred 2007 May 16 to 17 (MJD 54236.176-54237.813), 
lasting for 141 ks, with a total exposure for the XIS of 
114.8 ks after filtering. A 20.9 ks duration portion of this 
observation (MJD 54236.861-54236. 110) overlaps wit h 
the RXTE observations presented by IReig et all (2009) . 

The XIS observations were taken in 1 /4 window mode 
to avoid potential problems with pile-up. Our analysis 
began with the version 2.0.6.13 processing data (which 
has aspect correction for thermal wobbling) . The Charge 
Transfer Inefficiency (CTI) correction was remade and 
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Fig. 1. 
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light-curve of 4U 2206+54 from Suzaku XIS (XISO, XIS1 and XIS3 summed) in the 0.5 to 10 keV energy band. 



the event data cleaned using the standard filtering with 
the exception that the lower bound on the day side earth 
limb elevation was reduced from 20 degrees to 17.5 de- 
grees (which did not appear to increase the broad-band 
background). Events were selected from the active tele- 
scopes (XISO, XIS1, and XIS3) using 208" x 261" rect- 
angular extraction regions centered on the source, with 
the long axis aligned with the 1/4 window. As a limit 
on the background, events were also extracted from each 
telescope in regions of the same area as distant as possi- 
ble from the source. The average rate summed over the 
telescopes was 13.9 c s^ 1 for the source regions and 0.125 
cs^ 1 for the background regions. 

The XIS light-curve in the 0.5-10.0 keV range with 
~ 200 s binning is shown in Figure [TJ In order to avoid 
a loss of data, unequal bins sizes were used. Each good- 
time interval was divided into equal width bins, with the 
number of bins chosen to result in the bin width closest 
to 200 s. There appears to be pulsation with a period 
near 0.06 days, which is most evident in the spacing of 
the minima. Also evident is strong variablity at much 
longer timescales. 

2.1. Pulse Search Statistic 

To search for pulsation and determine its period we 
have developed a test statistic, which like the Lomb test 
statistic is based on hypothesis testing, but which can 
account for the low frequency variablity and other char- 
acteristics of the obse rvations. 

The Lomb method (|Lomblll975l ) , a standard technique 
for detection of pulsations in unevenly sampled data, as- 
sumes a null hypothesis of Gaussian white noise, and a 
alternate hypothesis of Gaussian white noise plus a sinu- 
siod. For 4U 2206+54 t he assumpt i on of Gaussian white 
noise is inappropriate. Rei g et all (|2009f) show that the 
power spectra of the flux in the 2-20 keV band are power- 
laws above 1 mHz, with an index between 1.5 and 1.9. 
The light curve in Figure [1] shows both short term flar- 
ing and strong long period (~ day) variability, suggesting 
this red noise spectra extends to lower frequencies. In ad- 
dition to the power-spectrum not being white, the distri- 
bution of the rates is far from being Gaussian. In the top 
panel of Figure [5] we show the distribution of the rates 
measured with Suzaku XIS, which is more exponentional 
than Gaussian. In the bottom panel we show the distri- 
bution of the natural log of the rate, which has a more 



symmetric distribution which is closer to a Gaussian. 

For a more realistic null hypothesis we will assume that 
the natural log of the flux is Gaussian noise, with a power 
spectrum of the form 
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where sq, a and T are adjustable constants and /o is a 
fixed reference frequency. The covariance of the natural 
logs of rates k and I is then 

poo 

V U = a 2 k S kl + / 5(/)sinc 2 (7r/r)cos(2^/[i fe -t ; ])4f (2) 
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where 07. is the error on = ln(rfc x 1 s), r is the average 

rate bin width, and the sine 2 term accounts for the effect 
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Fig. 2. — Distribution of the Suzaku XIS count rate (top panel), 
and of its natural log (bottom panel). 
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of binning. The likelihood of the null hypothesis is then 
given by 

Co = Det(27rF)-5 exp(-± J> fc - mKT^ - /"]) (3) 



where fj, is the mean natural log rate. In practice it is 
simpler for us to use the Cash statistic (jCash 1119791 ) C = 
— 21n£, giving 



Co = 5> fe - Hw - A*) - ln(Det(y)) 



ist. (4) 
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which is to be minimized with respect to sn, do, T and 

/x. 

For the alternative hypothesis the natural log of the 
flux is Gaussian noise with a power spectra of the same 
form, plus a sinusoidal profile. The resulting Cash statis- 
tic is 

Ci = Y,(yk-Pk)V kl 1 (m-pi)-ln{Vct(V)) + const. (5) 

kl 



where the pulse profile p, is given by 

Pk = M + a> c °s 27r/</j + b sin 2nftk 



(6) 



C\ is to be minimized with respect to sq, ao, T, /i, a, 
and 6 for a given trial frequency /. The test statistic 
analogous to the A% 2 is 



AC = min(Co) — min(Ci) 



(7) 



In the appendix we explain how the calculation of this 
statistic is implemented numerically. In Section 4 we 
describe how results from this statistic compare with the 
Lomb method. 

2.2. Suzaku Timing Results 

In Figure[3]we show AC for the Suzaku XIS light curve. 
The false detection probability of the highest peak is 5 x 
10~ 12 including the number of search trials within the 
0-1.0 mHz interval. 

We have investigated the presence of harmonics of the 
pulse frequency by including additional sinusoids in the 
profile (equation [6]) . Adding first harmonic terms de- 
creases the Cash C statistic by 9.9, which has chance 
probability of 0.7%. Adding second harmonic terms de- 
creases the Cash statistic an additional 11.4, which has a 

70 r 
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Fig. 3. — Delta Cash statistic versus frequency for the Suzaku 
XIS light curve. 
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Fig. 4. — 4U 2206+54 pulse profile obtained by epoch folding 
the Suzaku XIS 0.5-10.0 keV light curve with a period of 5554 s. 
The epoch of zero phase is MJD 54237.0448 (TDB). The smooth 
curve is the average of epoch folded simulated light curves based 
on the maximum likelihood fit. 

chance probability of 0.3%. Further harmonics continue 
to decrease the Cash statistic, but these decreases are 
not individually as significant. The period of 5554 ±9 s 
and power spectral parameters determined with these 
two harmonics included are listed in Table [TJ 

In Figure|4]we show the epoch folded Suzaku XIS rates 
using the period determined from the maximum likeli- 
hood fit with two harmonics included. The errors used 
are the square root of the sample variance in each phase 
bin divided by the number of rates in the bin. Using 
the pulse period, power spectrum, and profile parameters 
estimated from this same fit, we made 10000 simulated 
light curves with the same binning as in the Suzaku data. 
The simulation method is discussed in the appendix. The 
smooth curve in Figure[5]is the average epoch folded pro- 
file of these simulated light curves. 

In Figure [5] we show pulse profiles obtained from the 
Suzaku XIS and HXD/PIN data. Background sub- 
tracted light-curves with 24 s resolution were made for 
each energy band and epoch folded with a period of 
5554 s. For the PIN profiles light-curves where made 
from the cleaned event files, with live-time correction 
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Fig. 5.— 4U 2206+54 pulse profiles from the Suzaku XIS (on 
left) and HXD/PIN (on right). The epoch of zero phase is MJD 
54237.0448 (TDB). 
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TABLE 1 

Period and Power Spectra Parameter Estimates 



Observation 


Mid-Time 


Period 


So" 


a 


r 


Harmonics' 3 




(MJD) 


(s) 


(Hz" 1 ) 


(10~ 5 Hz) 






Suzaku XIS 


54237.0 


5554 ± 9 


31.6 ± 1.9 


2.2 ± 1.5° 


1.10 ±0.08 


3 


BcppoSAX MECS 


51141.1 


5420 ± 28 


21.4 ±3.6 


< 12 d 


1.14 ±0.18 


2 



a Power spectra normalization at frequency /n = lmHz. b Number of harmonics in fit including fundamental.' 
is 1.6 X 10~ 6 - 5.1 X 10" 5 Hz. d 90% confidence upper limit. 
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Fig. 6. — Power spectra estimated from a gapless 38 ks segment 
of the Suzaku XIS observations. The solid curve is average power- 
spectra of simulated light curves based on the maximum likelihood 
fit to the Suzaku XIS light curve. 



based on 144 s averages of the pseudo events. The sub- 
tracted background was from the version 2 non- X-ray 
background file plus the cosmic x-ray background calcu- 
lated with the flat response per the recipe on the Suzaku 
Data Analysis web page 6 . For both XIS and PIN pro- 
files the errors were calculated from the sample variance 
in each bin. Short flares in the data has resulted in the 
errors being correlated between neighboring bins. 

The modulation fraction of the profiles ([max- 
min]/[max+min]) is 49% in 0.5-2.5 energy band, increas- 
ing to 55% in the 20-30 keV band. A notch at the mini- 
mum is apparent in all profiles up to 30 keV. The profile 
is symmetric at low energies, with the leading edge be- 
coming stronger than the tailing edge at higher energies. 

We show in Figure [6] how the estimate of the pulse pro- 
file and power spectrum of the natural log of the rates 
is related to the empirically determined power spectrum 
of the rates at higher frequencies. For the higher fre- 
quency power spectrum we used a 37ks segment of the 
Suzaku XIS observations which contained no gaps. We 
made a light curve of this segment with 2 s resolution 
and constructed from it power spectra estimates based 
on the Fourier transform of the counts. We then pro- 
duced 10000 simulated lightcuves using the parameters 
estimated in the maximum likelihood fit of the Suzaku 

6 http:: /suzaku. gsfc.nasa.gov/docs/suzaku/aehp_data_analysis. html 



XIS data. Unlike the previous simulations, these light 
curves had uniform sampling covering the duration of 
the observation. The average power spectrum from these 
simulations is shown by the solid curve in Figure [6] 

3. ANALYSIS OF THE BEPPOSAX LIGHT-CURVE 

We ha ve reanalyzed t he Be ppoSAX observation pre- 
sented by |M_i_iE_et_in llll). The source was observed 
with BeppoSAX on 1998 November 11. We extracted 
events in the 1.65-10.0 keV range from the version 2 
merged MECS 2 and 3 files obtained from the ASDC 
Multi-Mission Interactive Archive, using a 4' radius ex- 
traction circle of centered on the source. As a limit on the 
background, events were also extracted from an annulus 
with the same area centered on a radius of 6'. The mean 
rate from the source region was 0.41 c s _1 while that 
of the background region was 0.014 c s" 1 , with much of 
this likely associated with the wings of the point-spread 
function. 

In Figure[7]we show the MECS light -curve with ~ 200 s 
binning of the 1.65-10 keV events. Again there appears 
to be pulsation with a period near 0.06 days, which is 
most evident in the spacing of the minima. 

We show a search for long-period pulsations in Figure 
[5] using the AC statistic with the same modeling used for 
the Suzaku observation. The highest peak has a proba- 
bility of occurring by chance of 5 x 10 -8 including the 
number search trials in the 0-1.0 mHz interval. From the 
location of the highest peak we estimate a pulse period 
of 5393 ± 28s. The peak near 0.36 mHz suggests a first 
harmonic. We investigated the presence of harmonics of 
the pulse frequency by including additional sinusoids in 
the profile (equation [6]) . Adding first harmonic terms 
to the pulse profile decreases the Cash statistic by 8.3 
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7.— light-curve of 4U 2206+54 from BeppoSAX MECS in 
65-10 keV range. 
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Fig. 8. — Delta Cash statistic versus frequency for the BeppoSAX 
MECS light curve. 

(1.6% chance probability). Improvements from adding 
additional harmonics were not as significant. A period 
estimate of 5420 ± 25 s is obtained from the fit including 
the first harmonic. This and the estimated power spectra 
parameters are shown in table Q] 

4. PERIOD DETERMINATIONS USING SIMULATED 
SUZAKU AND BEPPOSAX DATA 

We find small but sometimes significant differences be- 
tween the periods estimated with our maximum like- 
lihood analysis and those determined with the Lomb 
method. We also find differences between the period 
determined with the maximum likelihood method when 
harmonics are included with the fundamental in the fit 
and the period determined when only the fundamental 
is included. Here we present simulations which examine 
the extent to which these differences are expected. Our 
simulation method is presented in the appendix. 

In Figure [9] we present the results of a search for pul- 
sation in th e XIS light-curve using the Lomb method 
(|Lomblll975l ). For the Lomb power we are using twice the 
normalization given in iPress et al.l (fl992). A significant 
peak is evident near 0.18 mHz. Using the assumptions 
standard for the Lomb periodigram the chance probabil- 
ity of this peak to occur within the 0-1.0 mHz range due 
to noise is 2 x 10~ 31 . However these assumptions are 
not correct, as we have discussed, and this significance is 
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Fig. 9. — Pulse period search results from the Suzaku XIS data, 
showing the Lomb power versus the trial pulse period. 



over-estimated. This should be compared with Figure [3l 
Note the decrease of the power of the lowest frequency 
peaks relative to the main peak using the AC statistic. 

The pulse period estimated from the peak of the Lomb 
power is 5541 ± lis which is nearly consistent with our 
result of 5554±9s using the maximum likelihood fit with 
first and second harmonic included. We simulated 10000 
Suzaku light curves using the same times and binning as 
the light curve in Figure [TJ and the pulse period, power 
spectra and Fourier coefficients obtained from our maxi- 
mum likelihood fit for this light curve. For each we used 
the Lomb method to estimate the pulse period. The pe- 
riod distribution was consistent with a Gaussian with 
mean displaced by — 1.8 ± 2s from the simulated period 
of 5554 s, and a standard deviation of 17.6 ± 1 s. So 
for the Suzaku observations we conclude that the Lomb 
method is essentially unbiased, but with underestimated 
errors. The period measured with the maximum likeli- 
hood fit including only the fundamental was consistent 
within error to our final result from the fit including first 
and second harmonics. 

To perform the same investigation for BeppoSAX data, 
we made simulated light curves with the time binning 
of the MECS light curve using the estimated parame- 
ters from the maximum likelihood fit including the first 
harmonic. Using 10000 simulated light curves we found 
the distribution of the pulse period determined with 
the Lomb method was Gaussian with a mean displaced 
—23.8 ± 0.3 s from the simulated period (5420 s) and 
standard deviation of 32 s. The period obtained using 
the Lomb method using the actual data is displaced by 
—44 s from the period from our final maximum likelihood 
fit, consistent with the distribution observed in the sim- 
ulations. 

The period determined with the maximum likelhood 
fit of the BeppoSAX data using only the fundamental 
was displaced by -27 s from the period of 5420 ± 25 s de- 
termined from our final fit which included the first har- 
monic. Using 200 of the simulated light curves, we found 
the distribution of the periods determined with our maxi- 
mum likelihood method with only the pulse fundamental 
had a mean displaced from the simulated period (5420 s) 
by — 17 ± 2 s and a standard deviation of 28 ± 2 s. 

Using 200 simulated BeepoSAX light curves we found 
no bias in the periods estimated with our maximum like- 
lihood method with the fundamental and first harmonic 
included. However the standard deviation of the esti- 
mated periods was 28.4 ± 1.4 s which larger than the 
error (25 s) we obtained from the curvature of the Cash 
statistic versus frequency curve for the actual data. This 
is reflected in the error given in table [T] 

5. PERIOD DETERMINATIONS FROM THE EXOSAT 
OBSERVATIONS 

Observations of 4U 2206+54 where made with the 
Medium Energy (ME) proportional counters on the Eu- 
ropean Space Agency's X-ray Observatory EXOSAT. 
EXOSAT/ME observed 4U 2206+54 on three occa- 
tions; 1993 August 8, 1984 December 7, and 27 June 
1985. These observa t ions were originally presented by 
Saraswat & Apparao (1993). 

iReig et all (|2009( )analvsed the ME standard product 
light curves from these observations, which are available 
at HEASARC, and reported pulsation with a 5525 ± 30 
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Fig. 10.— The EXOSAT/ME light curves for the 1983 and 1985 
observations of 4U 2206+54 with 200 s binning. 

period. In figure ITD1 we show the light curves for the first 
and third observation with 200 s binning. We concluded 
that the background was not correctly subtracted for the 
standard products light curve of the second observation, 
and therefore have not shown it. 

Given the short durations of these light curves (9630 s, 
9400 s, and 7404 s respectively) relative to the pulse pe- 
riod, we do not believe a period measurement based on 
them is reliable. To demonstrate this we have made simu- 
lations based on two long gapless intervals of the Suzaku 
XIS light curve. For the first simulation we selected a 
number of overlapping 9600 s data segments within each 
gapless interval and from each segment made a pulse pe- 
riod e stimate using the Suzaku XIS rates in that segment. 
Since iReig et all (|2009h used the CLEAN method for pe- 
riod determination we did also. CLEAN attempts to it- 
eratively reconstruct the Four ier amplitude spectrum of 
the data (|Roberts et al.lll987l ). We used CLEAN within 
the program PERIOD version 5.0-2, distributed with the 
Starlink Software Collection. In our analysis of each seg- 
ment we used 5 iterations with a loop gain of 0.2. In 
the top panel of figure [TT] we show the period from each 
segment versus the midpoint time of the segment. The 
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11. — Periods determine using the CLEAN method using 
segments of the Suzaku XIS data. The top pannels are for 
segments and the bottom pannels for 7200 s segments. 



results of the second simulation, which used 7400 s seg- 
ments, is shown in the bottom panel. It is clear that 
the period estimates from these short data segments vary 
widely from the 5554 s period determined from the whole 
observations, and that the formal errors, which do not ac- 
count for the underlying red noise, grossly under estimate 
the error distribution. 

6. DISCUSSION 

We interpret the observed X-ray pulsations of 
4U 2206+54 as due to the rotation of an accreting magne- 
tized neutron star. An alternative interpretation, which 
was also discussed for the 10000 s pulsatio ns of the low 
lumin osity wind- fed system 2S 0114+650 (|Corbet et al.l 
1999), is that the observed period originates from a pul- 
sation of the companion BD +53° 2790. However, this 
explanation seems less plausible because such persistent 
modulations should be detectable in the optical photom- 
etry, but have not been reported. A mechanism for trans- 
ferring the modulations to the stellar wind that results 
in a large modulation of the x-ray flux has not been ad- 
vanced. In addition, the evolution of the pulse profile 
with energy that we see in the Suzaku data would not be 
expected. 

Is it reasonable for us to assume that such a long pe- 
riod accreting pulsar can be formed? The longest period 
which a neutron star can reach in its evolution is the 
period at which it em erges from the subsonic propeller 
state (|Ikhsanovll2507t ) 



Pbr = 15000 p 32 



16/21 -4/21 



-5/7 



10 15 gs" 



s, (8) 



where /132 and m are the dipole magnetic moment and 
mass of the neutron star in units of 10 32 Gcm 3 , and 
1.5 Mq, respectively. M c = ttRqpV 1c i is the mass with 
which a neutron star of a mass M ns interacts in a unit 
time moving through the wind of an average density 
p with a relative velocity V+ e i, and Rg — 2GM ns /V^ el 
is the Bondi radius. For the Suzaku observations we 
find a 0.5 — 70keV luminosity of 4 x 10 35 ergs" 1 im- 
plying a mass accretion rate onto the stellar surface of 
AI a = r DS L x/GM ns — 2 x lQ ^gs" 1 , assuming a distance 
of 2.6 kpc (Blay_eiaL|[2006). The mass accretion rates 
we inferred from published observations all fall in the 
range of 4 x 10 13 — 3 x 10 15 gs _1 , with the average near 
the higher end. Assuming M c ~ M a (i.e. direct accre- 
tion scenario) and putting this to Eq. [5J one finds that 
condition pbr <; 5500 s is satisfied only if the surface field 
of the neutron star in the present epoch is Bq > 10 14 G. 

An independent estimate of the field strength can be 
found considering the neutron star spin evolution. The 
average rate of the spin frequency change between the 
BeppoSAX and Suzaku observations was v = (— 1.7 ± 
0.3) x 10 _14 Hzs _1 . This indicates that the average 
spin-down torque applied to the star during this time 
was K s d > 2irl\i>\ : where / is the star's moment of 
inertia. If we now adopt the canonical prescription 
for the spin-down torque (K s d = ktn 2 /r l, see, e.g 



iLvnden-Bell k Pringld [1973 iLipunovl |l992|) , one finds 



A 1 ^ 



where 



10 32 k t 



-V2. 



1/2.1/2.1/2 / 
1 J 45 ^-14 I 



P 



5500 s 



Gem 3 . (9) 
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Here V- U = | z> | / 10 14 Hz s"\ and / 45 = 10 45 gcm 2 . 

1 /3 

r c = (GAfp 2 /47r 2 ) is the corotation radius of the neu- 
tron star, and fc t is a dimensionless parameter limited 
to fc t < 1. The value of \x m represents the lower limit 
to the dipole magnetic moment of the neutron star since 
the spin-up torque in the above calculations was assumed 
to be negligibly small. Therefore, our estimate remains 
valid independently of whether the star between the Bep- 
poSAX and Suzaku observations was persistently in the 
accretor state or its state was temporarily changed to the 
propeller. Thus, the 5500 s pulsations in the X-ray flux 
of 4U 2206+54 can be explained in terms of the spin pe- 
riod of the degenerate companion provided the neutron 
star in this system is a magnetar whose surface field at 
the present epoch exceeds 10 14 G. 

The modeling of 4U 2206+54 in terms of an accretion- 
powered magnetar suggests that the neutron star is rel- 
atively young (its age is limited to the characteristic 
time of the supercritical magnetic field decay), with an 
extended magnetosphere (r m cx /i 4 ^ 7 ), and a relatively 
small area of hot polar caps at the base of the accretion 
column (A p oc /x -8 / 7 ). However, the assumption about 
spherical geometry of the accretion flow, used for making 
the estimates ((8j) and ([9]) encounters major difficulties in 
explaining the mode by which the accretion flow enters 
the magnetic field of the neutron star at the magneto - 
sph eric boundary. A s first s hown bv lArons" fc Led ((19761) . 
and lElsner fe Lambl (| 19771) . a steady accretion onto the 
stellar surface in this case could occur only if the X-ray 
luminosity of the pulsar meets the condition L x > L cr , 
where 

6 ergs \ (10) 



io 37 ^ ™ 1/2 



Here re is the radius of the neutron star in units of 
10 6 cm. Otherwise, the interchange instabilities of the 
magnetospheric boundary are suppressed and the entry 
rate of accretion flow into the magnetosphere is too small 
to explain the X-ray luminosity of the pulsar. The star 
in this case wo uld appear rather as a burster than a per- 
sistent source (jLamb et al.lll977| ). 

It is easy to see, however, that 4U 2206+54 does not 
meet the above criterion. Nevertheless, it is a persistent 
accretion-powered X-ray pulsar with a relatively small 
amplitude of variations in X-ray intensity. This con- 
tradiction may indicate tha t either the spherical accre- 
tion is realized in a settling (|Naravan fc M edvcdey |2003t 
Ilkhsanovll2005l ) rather than direct accretion mode or the 
neutron star is accreting material from a disk. 

A scenario in which the neutron star accretes ma- 
terial from a hot envelope in a settling mode on the 
bremsstrahlung cooling t ime scale impli es the following 
condition to be satisfied (|Ikhsanovll2005l ) 



tbi{r m ) 
tg(r m ) 



GM n ' 



(11) 



Here iff(r m ) and ibr(?"m) are the free-fall time and 
bremsstrahlung cooling time at the magnetospheric 
boundary. Using parameters of 4U 2206+54 one finds 
M= ^ 10 17 gs _1 . Therefore, the persistent behavior of 
the source can be explained in terms of the settling accre- 
tion scenario only if the strength of the wind overflowing 
the neutron star exceeds the typical mass-transfer rate in 
a binary system with a 09.5V by almost four orders of 



magnitude (|Reig et al.ll2009D . Currently available obser- 
vations do not favor the assumption about so intensive 
outflow from the normal companion. However, if it were 
valid the condition pbr > 5500 s could be satisfied only 
for the surface field of the neutron star >4x 10 15 G. 

It is interesting that a similar value of the magnetic 
field can be found considering a situation in which the 
neutron star accretes material from a disk. A spin- 
down of the neutron star in this case indicates that its 
spin period is smaller than the equilibrium period, p eCl , 
which is defined by equating the acceleration, K su — 
A/VGAf ns r m , and deceleration, K sd — k t ^i 2 /r^, torques. 
Using parameters of 4U 2206+54, one finds that the con- 
dition p < p oq is satisfied only if the strength of the dipole 
field at the stellar surface in the present epoch is limited 
to 



B(r ns ) > 3x 10 15 nlTK 7/12 m 5 / 6 MUl 



V 5500s/ 



7/6 



G, 
(12) 

where Ko.5 = re/0. 5 is the parameter accounting the 
geometry of the acc retion flow normalized following 
iGhosh fc Lanlbl (fl97l . The spin-down time scale of a 
newly formed neutron star to a period of 55 00s under 
these conditions is close to 5000 yr (see Eq. 20 in llkhsanovl 
|2007| ). This exceeds the decay time of the m agnetic field 
confi ned to the crust of the neutron star (jColpi et al.1 
2000), but is still smaller than the decay time of the 
field permeating the core. Since the X-ray spectrum of 
the source resembles spectra of accretion-powered pul- 
sars one can assume that the contribution of the source 
associated with the field decay to the system X-ray lu- 
minosity is small. The field decay time in this case can 
be limited to 



i doc > 2 x 10 5 ^3.3 r 6 3 L 35.6 y r > 



(13) 



which is consistent within results of iHevl fc Kulkarnil 
(1998), who treated the magneto-thermal evolution of 
the star in terms of ambipolar diffusion. Here L35. 6 
is the X-ray luminosity of the source express in unites 
of 4 x lQ 35 ergs _1 . The proton cyclotron feature in 
this case is expected to be observed at the energy of 
~ 30 (B /5 x 10 15 G) keV. 

It, therefore, appears that our interpretation of the 
observed 5500 s pulsations in terms of spin period of the 
neutron star is reasonable and suggests that the neu- 
tron star in 4U 2206+54 is a magnetar with a surface 
field of about (3 — 5) x 10 15 G. A question about the 
geometry of the accretion flow remains, however, open 
so far. The spherical approximation of the flow geome- 
try can be applied only under the condition M c 3> M a . 
On the other hand, the observed wind velocity of the O - 
star companion, V w > 350 kins" 1 , (see lRibo et al.f 2006'). 
substantially exceeds the upper limit to the relative ve- 
locity between the star and the wind at which the disk 
formation could be expected, 

V rcl < V CI ~ 120^/ 2 4 Ma/ 714 m"/ 28 MT P2o /4 kms- 1 . 

(14) 

Here P20 is the orbital period of the system in units of 
20 days, and £0.2 = £/0.2 is the parameter accounting 
for the inhomogeneit ies of the accr etion flow, which is 
normalized following iRuffer t (1999). Nevertheless, one 
cannot exclude a presence of a fossil disk surrounding the 
magnetosphere of the neutron star which could be formed 



Finger, Ikhsanov, Wilson-Hodge & Patel 



at the fall-back stage of the magnet ar evolutionary track 
(|Wooslevlll98i lEksi k Alparll2003f ). Further studies of 
appearance of such a disk in a situation when the neutron 
star rotates extremely slowly and has a huge magneto- 
sphere could help to choose appropriate flow geometry 
and reconstruct the accretion picture in this enigmatic 
source. 

Finally, we would like to point out here that at the 
rate v = (-1.7 ± 0.3) x lCT^Hzs- 1 the pulsar would 
reach zero frequency in 300 yr, making it is unlikely that 
this is a long-term trend. Over long time-scales wind- 
fed pulsars tend to show random-walk behavior in their 
spin frequency. For 4U 2206+54 we can estimate the 
random walk strength as S = < {Av) 2 / At > ~ 8 x 
10 _20 Hz 2 s _1 . For compa rison Vela X-l has a random 
walk strength of 2 x 10^ 20 (jDeeter et al.lll987D . The ran- 
dom walk behavior must break down, and the torque 
become correlated at short time-scales since the spin-up 
rate v ~ (S/At) 1 ^ 2 from wind accretion is limited by 
that due to disk accretion at the same mass accretion 



rate. Assuming a mean mass accretion rate of 10 15 g s" 1 , 
and a magnetic field of 10 15 G, we find the spin-up rate 
must be correlated on timescales shorter than ~ 5 days. 

Based on the s trong emi ssion in Helium lines 
iNegueruela fc Reid ()2001l ) and iBlav et all (|2006) have 
proposed that the optical companion BD +53° 2790 is 
similar to 9 1 OriC, which has a oblique magnetic dipole 
with kG surface strength which magnetically channels its 
stellar wind. If so, this long correlation timescale could 
be related to long-term correlations in the wind speed or 
density associated with this magnetic channeling. 
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APPENDIX 

CALCULATION OF THE CASH STATISTICS C AND d 
Numerically we calculate V in equation [5] as 

M 

Vu = a 2 k 5 M + J2 w J -5(/ i )siiic 2 (7r/ i r) cos(27r/ J -[t fc - U])Af . (Al) 

3=0 

Here the frequency fj = jAf and the frequency step used is A/ = (16T) -1 with T the duration of the observation. 
The weigths wj give trapazoidal integration with wj — 1 except for wq — wn = 0.5. M is chosen as the nearest integer 
to T/t so that the integral extends to where the power is cut off by the binning effect. 
By using 

cos(27r/j [tk — ti)) = cos 2nfjtk cos 2irfjti + sin 2itfjtk sin 2itfjti (A2) 
equation IA1I can be express as 

V = U T U (A3) 

where the matrix U has N + 2M + 1 rows and N columns, with N is the number of measurements, and T signifies 
the matrix transpose. Explicitly U is given by 



Ui k = a i 5 i,k for < i < N 

U N+2 j-x,k = K-5(/j)sinc 2 (7r/j-r)A/] 2 cos2w fjt k for l<j<M 
U N +2j,k = [u; 3 5(/ J )sinc 2 (7r/ J T)A/]^sin27r/^ fe for 1<3<M 



By applying a series of N orthgonal Householder transformations (reflections) U can be transformed into an N x N 
upper triangular matrix U, while perserving equation I A3I (|Householderlll958D . We then have 

Ci = |E/- T (y-p)| 2 +2^1n(|^|), (A5) 

i 

with a simular expression for Cq. Here y is the vector with components yk and p is the vector with components pk- 
We will introduce the augmented parameter vector q = [fi, a, 6,-1], and using equation [5] define the matrix H so that 
p y = Hq. Then with R = U~ T H we have 

d = |i?q| 2 +2j2H\Uu\) , (A6) 

i 

By applying Householder transformations R can be transformed into a 4 x 4 upper triangular matrix R while preserving 
equation I A6I We can then minimize C\ with respect to /x, a and b by setting the first three components of i?q to zero 
and solving for /i, a and 6, with the minimum of |i?q| 2 given by R 2 3 . 
We minimize Co and C± with respect to the power spectral parameters sq, qq, and Y by the downhill simplex method. 
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LIGHT CURVE SIMULATIONS 

In our simulations we first simulate natural log rates with covariance given by equation l All with the measurement 
errors a k set to zero, and means pk given by equation [6] (perhaps with terms for harmonics). Simulated counting 
statistical errors are then added to the rates r k - 

As above an upper triangular matrix U is calculated such that U T U = V. Then with x is vector of N random 
normal varients with unit variance zero mean we calculate 

* = U T x (BI) 

The covariance matrix of '3/ is 

< >= U T < xx T > U = U T U = V . (B2) 

The simulated rates prior (without counting statistical errors) are given by 

r fc = exp(tf fc +p fc ) . (B3) 
Counting statistical errors are then simulated from these rates and the bin widths. 
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